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Abstract 

In a recent paper we proposed the expansion of the space of variations in energy calculations by 
considering the approximate wave function ip io he a functional of functions X • = rather 
than a function. For the determination of such a wave function functional, a constrained search is 
first performed over the subspace of all functions x such that ■(/'[x] satisfies a physical constraint or 
leads to the known value of an observable. A rigorous upper bound to the energy is then obtained 
by application of the variational principle. To demonstrate the advantages of the expansion of 
variational space, we apply the constrained-search-variational method to the ground state of the 
negative ion of atomic Hydrogen, the Helium atom, and its isoelectronic sequence. The method 
is equally applicable to excited states, and its extension to such states in conjunction with the 
theorem of Theophilou is also described. 
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1. Introduction 

In the traditional application of the variational principle for the energy Q|, the space of 
variations is limited by the choice of the analytical form chosen for the approximate wave 
function. For example, if Gaussian or Slater-type orbitals or a linear combination of such 
orbitals is employed in the energy functional of the wave functions, then the variational 
space is limited by this choice of functions of the wave functions. In a recent paper ^2] we 
proposed the idea of overcoming this limitation by expanding the space over which the 
variations are performed. This allows for a greater flexibility for the structure of the wave 
function. A consequence of this greater variational freedom is that better energies can be 
obtained. Or, equivalently, fewer variational parameters are needed to achieve a required 
accuracy. 

The manner by which the space of variations can be expanded in principle is by 
considering the approximate wave function ^/l) to be a functional of a set of functions 
X : ^/^ = rather than a function. The space of variations is expanded because the 

functional ip[x] can be adjusted through the function x to reproduce any well behaved 
function. The space over which the search for the functions x is to be performed, however, 
is too large for practical purposes, and a subset of this space must be considered. The 
subspace over which the search for the functions x is to be performed is defined by the 
requirement that the wave function functional ip[x\ satisfy a constraint. Typical constraints 
on the functional ipix] ^''^^ those of normalization, the satisfaction of the Fermi-Coulomb 
hole sum rule, the requirement that it lead to observables such as the electron density, 
diamagnetic susceptibility, nuclear magnetic constant, Fermi contact term, or any other 
physical property of interest. With the wave function functional ip[x] thus determined, a 
rigorous upper bound to the energy is obtained by application of the variational principle. 
In this manner, a particular property of interest is obtained exactly while simultaneously 
the energy is determined accurately. We refer to this way of determining an approximate 
wave function as the constrained-search-variational method. The method is general in that 
it is applicable to both ground and excited states. 

In section 2 of the paper, we explain the constrained-search-variational method in 
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further detail. To demonstrate the ideas involved, we apply the method in section 3 to the 
ground state of the Helium atom, its isoelectronic sequence, and the negative ion of atomic 
Hydrogen. Concluding remarks are made in section 4. 



2. Constrained- search variational method 

To explain the method for the determination of a wave function functional, consider the 
non-relativistic Hamiltonian of the Helium atom, the ions of its isoelectronic sequence, and 
the negative ion of atomic Hydrogen. In atomic units {e = h = m = 1) 

^ 1 o 1 o Z Z \ . . 

H = -:tV - -V + — , 1 

2 2 ri r2 ri2 

where ri, r2 are the coordinates of the two electrons, ru is the distance between them, 

n 

and Z is the atomic number. In terms of the Hylleraas coordinates ^Si]: s = ri + r2, t = 
ri — r2, u = ri2, which are the natural coordinates for this atom, we choose the approximate 
wave function functional to be of the general form 

ij[x] = Hs,t,u)[l-f{x]s,t,u)], (2) 

with $(s, t, u) a Slater determinantal pre-factor and /(x; s, t, u)] a correlated correction term: 

/(s, t, u) = e-'^"(l + qu) [1 - s, t, u){l + n/2)], (3) 

where g is a variational parameter. Note any two electron wave function in a ground or 
excited state maybe expressed in this form. The Slater determinant may be chosen to be the 
Hartree-Fock theory wave function jj], or determined self-consistently within the framework 
of Quantal Density Functional Theory |5]. For purposes of explanation, we consider here 
the determinant composed of Hydrogenic functions. Thus, for the ground state l^S* of 
the Helium atom we have = (a^/7r)e~°'', and for the excited triplet 2^5* state 

s,t] = ^J2J^{a^ I ■n)e~^^t. (In the latter, for explanatory purposes, screening effects are 
ignored). Further, we assume that x is a function only of the variable s: \E' = s)]. 



The approximate wave function 
electron-electron cusp condition 



unctional \E'[x(q', s)] for the ground state then satisfies the 
It also satisfies the electron-nucleus cusp condition for 
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a = Z. 



Next consider observables such as the size of the atom, diamagnetic susceptibihty, nuclear 
magnetic constant, Fermi contact term, etc, which are represented by the expectation of 
operators W = ri + r2, W = rf + r^, W = 1/ri + l/r2, W = 5(ri) + 5(r2), respectively. For 
the normalization constraint W = 1. In terms of the Hylleraas coordinates, these operators 
are W{s) = s, W{s,t) = {s^ + t^)/2, W{s,t) = W{s,t) = ^[^J^ + ^J^], and 

W{s, t) = 1. In general, observables can be represented by single-particle operators expressed 
as W{s,t). The expectation of the operator W{s,t) is then 

(W) = j ^*[x]W^x]dr = (W^o) + Aiy, (4) 

where (for the ground state) 

{W)^ = j |$(a,s)|2l^(s,t)dr, (5) 

AW = j ma,s)\^W{s,t)[f\q,s,t,u)-2f{q,s,t,u)]dT (6) 

POD 

= ma,s)\'gis)ds, (7) 



^0 

where 

g{s)= udu dtW{s,t){s''-e)[f{q,s,t,u)-2f{q,s,t,u)]. (8) 

JQ Jo 

We now assume that the expectation (W) is known either through experiment or via 
some accurate calculation 1?^. As our choice of $(a;; s) is analytical, then both {W)o and 
AW are now known. 



The next step is the constrained search over functions x{Qi for which the expectation 
(W) of Eq.(4) is obtained. If the parameter a in Eq.(7) is fixed, then there exist many 
functions g{s) for which the expectation (W) can be obtained. This corresponds to a large 
subspace of wave function functionals (See Ref. 2). On the other hand, if the parameter a 
is variable, then the only way in which Eq.(7) can be satisfied is if 

g{s) = G, (9) 
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where G is a determinable constant. This is equivalent to the constrained search of all wave 
function functionals over the subspace in which Eq. ( 7) is satisfied. 

As an example consider the normalization constraint for which (W) = {W)q = 1, so that 
AW = 0. Then the only way in which Eq.(7) can be satisfied, (for variable a) is if 

g{s) = 0. (10) 

This condition is thus equivalent to the constrained search over the subspace of all 
normalized functionals ^'[x(g, s)]. 

Substitution of f{x',s,t,u) into Eq.(lO) leads to a quadratic equation for the function 
X(?,s): 

a((l, s)x{(l, sf + ^biq, s)xiq, s) + c(q, s) = 0, (11) 

where 

a{q, s) = [\s^u^ - uy3){l + ^2)^(1 + qufe-^'^'^du, (12) 
Jo 

%,5) = - / (sV-M73)(l+M/2)(l + gM) 

Jo 

[e-^i^{l + qu)-e-'i'']du, (13) 

c(g,s)= / (sV-KV3)(l + g^i)[e"'''"(l + gH)-2e-«"]dii. (14) 
Jo 

The integrals for the coefficents a{q,s), b{q,s), and c{q,s) are determined analytically. 
Solution of the quadratic equation is equivalent to searching over the entire subspace of 
normalized wave function functionals. In this example, the subspace corresponds to two 
points. The two solutions Xi(g; s) and X2(?; s) lead to two normalized wave functions ip\xi] 
and V'[X2]. 

The generalization to the case when W = W{s) or W = W{s,t) follows readily. In 
either case, one has also to solve a quadratic equation for the determination of the functions 
x(?, a; s). One thus obtains two wave function functionals that lead to the exact value for 
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{W). 



For the normahzed wave function functionals determined above, the energy functional in 
terms of (s, t, u) coordinates which is 

ilj*HiljdT (15) 



= 27,' j^ds fdu fdtM^"-t")[i^y+i%y+ 

-[AZsu - is' - f)]^p'}, (16) 

is then minimized with respect to the parameters a and q. (The prefactor minimizes the 
energy at a = Z — 5/16). 

For wave function functionals determined by sum rules other than normalization, the 
functional /['^/'[x]] must be divided by the normalization integral J ip*ipdT. In this manner, 
the wave function functionals '4'[x\ ^^e normalized , obtain the exact value of the expectation 
{W{s,t)), and lead to an accurate value for the ground state energy. 



3. Application to two-electron atomic and ionic systems 

In this section, we apply the constrained-search-variational method to two-electron 
atomic and ionic systems. The two wave function functionals ip[xi\ and ip[X2\ employed are 
those determined via the constraint of normalization as described in the previous section 
with the crude Hydrogenic prefactor. In Table I we quote the values for the ground state 
energy for H~ , the He atom, and its isoelectronic sequence. For the He atom we also quote 
;he values of Hartree-Fock theory jj], the 3-parameter Caratzoulas-Knowles wave function 
^, and the 1078-parameter Pekeris wave function 0. For H and the other negative ions 
corresponding; to Z = 3 — 8, we give the values of the variational-perturbation results of 
Aashamar j9|. The satisfaction of the virial theorem and the percent errors as compared 
to the Pekeris and Aashamar values are also given. The functions Xi{Qi s) and X2{(li s) for 
B^+, and 0^+ are plotted in Figs. 1-3. 
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FIG. 1: Fig.l: The functions Xi{Q:s) and X2{q,s) for H 



Observe that the improvement of the energies of the two wave function functional over 
the prefactor values is generally an order of magnitude. As expected, the energies as well 
as the satisfaction of the virial theorem improves with increasing atomic number Z. For 
the He atom, the energies of both V'fXi] and ip[x2] are superior to those of Hartree-Fock 
theory and of the 3-parameter Caratzoulas-Knowles wave function. Furthermore, whereas 
the prefactor leads to a negative electron affinity, both wave function functionals lead 
to a positive electron affinity for H~ as must be the case as H~ is stable. The exact 
satisfaction of the virial theorem by the prefactor is a consequence of scaling, whereas that 
of Hartree-Fock theory is because of self-consistency. 

In Table II we quote the values of the operators W — = —2,-1,1,2, and 
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TABLE I: Rigorous upper bounds to the ground state of R- , He, Li+, Be^+ , B^+, C^^, N^+, 
O^^, in atomic units as determined from the wave function functionals determined via the con- 
straint of normahzation together with the values due to Hartree-Fock (HF) theory |^], Caratzoulas- 
Knowles(CK)y|, Pekerisj?] and Aashamar Q]. The satisfaction of the virial theorem, and the per- 
cent errors compared to the values of Pekeris and Aashamar are also given. 



Ion or Atom 


Wave function 


Parameters 


Ground state energy 


% error 
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2.0019 
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-0.52775 
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Xi 
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—22 01729 
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X . C X 




X2 




a = 4.6592,9 = 0.4401 


-22.01973 


0.051 


1.9997 


Aashamar 




-22.03097 




2.0000 








-32.34766 


0.181 


2.0000 




Xi 




a = 5.6574, g = 3.2528 


-32.39230 


0.043 


1.9993 




X2 




a = 5.6578, g = 0.4839 


-32.39511 


0.034 


1.9997 


Aashamar 




-32.40625 




2.0000 
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-44.72266 
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2.0000 
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Xi 




a = 6.6572, q = 3.9381 


-44.76729 


0.032 


1.9995 




x.2 




a = 6.6584, g = 0.5511 


-44.77035 


0.025 
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Aashamar 




-44.78145 




2.0000 
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1.9998 


Aashamar 




-59.15660 




2.0000 
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FIG. 2: Fig. 2: The functions Xiil^^) and X2{q,s) for i?^"*". 



ly = S{ri) + 5(r2) for the He atom as determined by both ip[xi\ and 'ip[x2] together with 
the Hartree-Fock theory, Caratzoulas-Knowles, and Pekeris values. The accuracy of these 
results is, of course, not correct to second order as are those for the energy. Nonetheless, the 
results are considerable improvements over the prefactor values. They are also all superior 
to the 3-parameter results of Caratzoulas-Knowles. The latter indicates that the two 
wave function functionals although also determined via energy minimization, are superior 
throughout space. Thus, by expanding the space of variations, one obtains a superior wave 
function not only in the region contributing most to the energy, but also in other regions of 
space. The superiority of the Hartree-Fock theory values, on the other hand, is due to the 
fact that in this theory, the expectations of single-particle operators is correct to second 
order jl^ . 
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Wc note that the two functions Xi(g, s) and X2{(l,s) are very different from each other: 
Xi(g, s) is positive whereas X2{(1, s) is negative. See Figs. 1-3. Thus, although the analytical 
form of the wave function is the same, the two functional ■0[xi] '0[X2] are very 
different. Nevertheless, they lead to accurate results that are essentially the same. Thus, 
the constrained search for the functions x over this subspace of normalized wave functions 
leads to two physically meaningful functionals. 



4. Concluding remarks 

In this paper we have shown how to expand the space of variations in calculations of 
the energy by constructing approximate wave functions that are functionals rather than 
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TABLE II: The expectation value of the operator W = Ya=i n = -2, -1, 1, 2 and = 5{ti) + 
(5(r2) for the He atom employing the wave function functionals determined by the normalization 
constraint, and by the Hartree-Fock theory(HF)0], Caratzoulas-Knowles(CK)y|, and Pekeris [3| 
wave functions (WF). 



WF 


< (1/ri + l/r2) > 


< (l/r? + 1/ri) > 


< {rj + rl) > 


< ri + r2 > 


Wri)+'^(r2)) 


^> 


3.3750 


11.391 


2.1069 


1.7778 


3.05922 


V'[xi] 


3.3773 


11.726 


2.1924 


1.8057 


3.37921 




3.3784 


11.727 


2.1876 


1.8041 


3.37925 


HF 


3.3746 


11.991 


2.3697 


1.8545 


3.5964 


CK 


3.3911 


11.714 


2.1292 


1.7848 




Pekeris 


3.3766 


12.035 


2.3870 


1.8589 


3.6208 



functions. The wave function functionals depend upon functions that are chosen so as to 
satisfy a sum rule or reproduce the value of an observable. In this constrained-search- 
variational method, wave functions that are accurate over all space are thereby obtained. 
The framework presented is general and applicable to both ground and excited states. For 
excited states, one would in addition employ the theorem of Theophilou according to 
which if ipi,ip2, ,fm,---y cire orthonormal trial functions for the m lowest eigenstates of the 
Hamiltonian H, having exact eigenvalues Ei, E2, Em,--- , then Yl^iiy^il^lfi) — ■ 
In this way, a rigorous upper bound to the sum of the ground and excited states is achieved. 
With the ground state energy known, a rigorous upper bound to the excited state energy is 
then determined, while simultaneously a physical constraint or sum rule is satisfied or an 
observable obtained exactly. 



In the calculations presented to demonstrate these ideas, a crude Hydrogenic Slater 
determinantal prefactor was employed. Improved results may be obtained through a better 
prefactor. Fully self-consistently determined prefactors for many-electron systems may 
be achieved, for example, via Quantal density functional theory. The latter is a local 
effective potential energy theory of noninteracting Fermions with the true density in which 
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the multiplicative potential energy operator representative of all the many-body effects 
is explicitly defined in terms of the interacting system wave function and the orbitals of 
this model system. These orbitals, determined self-consistently, then constitute the Slater 
determinantal prefactor. Or one could employ analytical or self-consistently determined 
Hartree-Fock theory orbitals for the prefactor. Another step towards improved results 
would be to further expand the space of variations defining the functions x- such a 
case, the equation for the functions x could be an integral equation. Other analytical 
forms for the correlation factor could also be employed. These avenues are being pursued 
to better understand the ideas underlying the construction of wave function function- 
als, and to employ these functionals within the context of Quantal density functional theory. 
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